NN Annihilation in Large iV c QCD with p and uj> Mesons 



Yang Lu[| and R. D. AmadoQ 
Department of Physics, 
University of Pennsylvania, Philadelphia, PA 19104, USA 

(May 16, 1995) 
Abstract 

We use classical/large iV c (number of colors) QCD to study nucleon-antinucleon anni- 
hilation at rest. We include pion, rho and omega fields in the classical dynamics, starting 
from the nonlinear sigma model. We begin with a spherical blob of pionic matter with 
zero baryon number and energy of twice the nucleon mass. This evolves according to the 
classical dynamics to pion and vector meson fields into the radiation zone. The radiation 
fields are quantized using the method of coherent states modified to include isospin and 
four momentum conservation. Empirical information is extracted from that coherent state. 
We find good agreement with data for single particle momentum and number spectra, and 
with branching ratios to many annihilation channels. All this emerges with nearly no free 
parameters. 
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1 Introduction 



There is growing interest in using classical versions of QCD to model the dynamics in strongly 
interacting processes leading to many pions. The notion is that the many pion state can be 
approximated as classical pion radiation. Classical QCD is QCD in the limit of a large number 
of colors |l|, ||], (Nq), and is naturally expressed in terms of the classical pion field. The strong 
interaction dynamics is much easier in the classical than in the quantum guise, and the pion 
radiation can be quantized in the radiation zone using the method of coherent states. This 
general approach of classical dynamics leading to pion radiation that can then be quantized 
has been applied to the evolution and decay of the disoriented chiral condensate ||, to pion 
production in heavy ion collisions and to nucleon antinucleon annihilation 

Though no one has actually derived a classical pion field theory directly form QCD, it is gen- 
erally agreed that any such theory will feature baryons as nonperturbative topological "knots" 
of the pion field and corresponding topologically stable baryon number. Skyrme |7]] developed 
just such a theory long before the advent of QCD. Some of its details may be special to Skyrme's 
theory, but its general features should pertain to any classical manifestation of QCD. It is this 
Skyrme approach that we have taken to describe low energy proton antiproton annihilation. 
We have shown that the pion momentum, number and charge spectra emerge correctly from 
this picture § and that extending the Skyrme treatment to include the classical omega field 
makes it possible to also account for the omega mesons seen in low energy annihilation 

In this paper we extend the large Nc picture further to include the rho meson and study its 
effects on annihilation. We concentrate on the observables of nucleon antinucleon annihilation 
at rest and show that with very few parameters we are able to account for the major features 
of the pion momentum and number spectrum, the branching ratios into the many annihilation 
channels and the ratio of direct pions to those from vector meson decay H. Two principal 
features of our treatment seem to be crucial to giving this agreement. First, classical QCD 
generates the correct fraction of vector mesons during the annihilation process. We do not 
have to assume the vector mesons were present in the baryons prior to annihilation. Second, 
the coherent state, with the constraint of isospin and four momentum imposed, describes all 
the annihilation channels in terms of a single quantum state, with each channel obtained from 
projections on that state. Thus the different channels have a common origin and the branching 
ratios come out naturally from that state with the constraints of phase space already imposed. 

Early studies of annihilation in the Skyrme approach showed that when a Skyrmion 

and antiSkyrmion annihilate the classical pion field streams away from the annihilation region 
as fast as causality will permit, leaving in a coherent burst of energy. Thus far no one has 
used the Skyrme approach to model the colliding of a proton and antiproton and its subsequent 
radiation. There are no real conceptual obstacles to doing such a calculation but it is technically 
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challenging, and it has not been done. It is, of course, still far easier than the corresponding 
quantum calculation in lattice QCD. Our models of annihilation have tried to exploit the rapid 
appearance of the pion pulse and have studied the development of the pion wave after the 
annihilation region has been formed. We begin with a "blob" of pionic matter with total 
energy of two nucleons, zero baryon number and in a region of about one Fermi radius. For 
simplicity we take this region to be spherically symmetric. We then use the large Nq dynamics 
to propagate the pion field outward into the radiation zone, where the nonlinearities of the 
dynamics can be neglected, and use the classical pion wave in the radiation zone to construct a 
coherent state. The observed features of the pions from annihilation are then studied in terms 
of that state. To obtain agreement with experiment it is necessary to impose both isospin and 
energy-momentum conservation on the coherent state. We have done all this both in the pure 
Skyrme picture, which includes pions only ||, and in a generalized form that includes omega 
mesons ||. 

In this work we extend the theory to include rho mesons. They are included in the large 
Nq dynamics in a standard way [11]. We also focus on the phenomenology of annihilation. 
We find agreement with the major trends of branching ratios and spectra with essentially no 
free parameters. Thus including the rho meson not only yields considerable phenomenological 
improvement but it also widens the circle of QCD dynamics that can be modeled in this clas- 
sical/coherent state way. In the next section we briefly outline our plan of attack. Section 3 
gives the massive Yang-Mills formulation for including rho and omega fields in the Skyrme or 
large Nc dynamics. Section 4 details the initial state assumptions we make, how the asymptotic 
fields are obtained and their form. Section 5 gives the quantum coherent state calculation based 
on those asymptotic classical fields. Section 6 gives our results and Section 7 presents some 
conclusions and prospectives for future work. 



2 Plan of Attack 

In low energy pp — > rnr about 33% of pions come out not directly from the annihilation but 
from the decay of meson resonances fl2|| . These are mostly the vector meson resonances, p and 
lo. To introduce the vector mesons in large N c QCD, we treat them as massive Yang-Mills fields 
which gauge the Uy (2) symmetry of the non-linear a model (contains tt only). This non-linear a 
model is the starting point of the Skyrme model or of any other classical or large Nc treatment 
of QCD. The chiral anomaly gives the coupling between ujpir, responsible for the interactions 
to —* 7T + 7r - 7r°, lo — > pir etc. The SU(2) gauge coupling leads to p — > 2tt. Vector dominance is in 
the model and the baryon is stabilized by the vector meson terms, eliminating the need for the 
Skyrme term. Though no one has derived classical QCD from quantum QCD, it is generally 
believed that the classical picture will begin with the non-linear sigma term, include the gauged 
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vector mesons and continue in a series of terms with higher and higher derivatives. These 
higher terms become important at high momentum or short distance, but since annihilation at 
rest is, in some sense, a low energy process, we believe keeping only the first few terms in that 
expansion is a reasonable starting point. We further believe our results give some support to 
that view. 

With our classical field theory, we study the dynamics of annihilation at rest by starting 
with a "blob" of matter of baryon number zero, of size about 1 fm and of energy twice the 
nucleon mass, made up purely of pions. We use the dynamical equations to study the space- 
time evolution of this blob including the generation of classical vector meson fields. At large 
t and r, the fields decouple and we have free field solutions. We quantize these free waves by 
coherent state techniques. We project the coherent state onto states of good four-momentum 
and good isospin. 

This quantum coherent state, generated from our dynamical, classical field solutions, is the 
starting point of our phenomenology. We obtain the particle spectra and branching ratios from 
the coherent state simply by taking the appropriate matrix element projection of that state. 
This means that all the various annihilation channels originate from that single state and are 
thus related. This relationship coupled with the correct phase space imposed by four-momentum 
conservation, leads naturally to a unified and satisfactory account of the data. 



3 Massive Yang-Mills Fields: the Classical Solution 



We introduce the p and uj classical fields as C/(2)y-gauge fields in the non-linear a model. The 
Lagrangian for the non-linear a model is 



f 2 



(1) 



where U is a unitary field in isospin space and r space and where f n is the pion decay constant. 
This non-linear a term is also the leading part of the Skyrme Lagrangian @ . We introduce the 
p and u fields as massive gauge fields through the U(2)y covariant derivative 



(2) 



where r is the Pauli matrices for isospin and greek indices specify space-time components. 
(Note that isovectors are boldface while arrows and "hat" denote spatial vectors.) Since the u> 
field commutes with U, it does not appear in Its couplings are given entirely by the gauged 
Wess-Zumino anomaly term. The complete Lagrangian then is [11, O] 
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where £ CT is the gauged nonlinear a model with a term for the pion mass. The field strength 
tensors of the vector mesons are given by 
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The mass term for the vector mesons is C m and we relate the "charge" g to the pion decay 
constant and the vector meson mass (assumed to be the same for p and to), m = 770 MeV by 
the KSFR relation [14], m = The gauged Wess-Zumino term gives Twz- The baryon 

current is given by 



1 



24vr 2 



n»ap Tr \(u*d v U)(tfd a U)(UidpUj 



(8) 



and the Wess Zumino term for the pion field Y^ Z {U) = for S'L r (2)-valued U field. 

We now wish to build in our initial condition of a spherically symmetric blob of pionic 
matter. That means that the pion field will always be of the hedgehog type, and we may write 



U = exp [it • fF(r, t)] 



0) 



where F is essentially the pion field. Furthermore we can decompose the vector meson fields 
into a time and space components as follows 



w (r,t) p = 
f % u) r (r, t) p = i 



and p 



pa 



f a H(r, t) p = 

gr r° 



(10) 



Here a is an isopin label and p a space-time label. The spatial unit vector is denoted by f 
and the indices a,i,j,k range over 1,2,3. We have explicitly removed a 1/r from the spatial 
part of the rho field. With our initial conditions, the H field decouples and we will drop it in 
subsequent equations. In terms of these quantities, we can write 
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\ dr dt J 

The nonlinearities of the theory and the nature of the field couplings are clear in the various 
terms. It should be stressed again that though we have put c the speed of light to one in these 
equations, we have not had to do that for fi. These are purely classical equations, and there is 
no h. 

The equations of motion that follow from this Lagrangian are 
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The corresponding energy density is: 
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We use these equations to propagate our initial field configuration into the radiation zone. 
In the next section we turn to a discussion of that configuration. 



4 Initial Field Configurations and Asymptotic Amplitudes 

As in ||, we model the initial pion field as 

F(r, t = 0) = h T exp(-r/a) (22) 

r z _|_ a z 

with a a range parameter and h adjusted so that the total energy equals to twice the nucleon 
mass. In our previous studies we arbitrarily fixed the range by a = l/m n . Here, since we are 
interested in a better treatment of the phenomenology, we take a as a parameter and fit it to 
obtain the correct asymptotic one pion momentum distribution. We find that a = 1.1 fm leads 
to an excellent account of the observed pion momentum distribution, (see Section 6) The initial 
field configurations are assumed to be static, so that at t = 

F = LO r = dj = G = 0. (23) 

We also take the initial vector meson fields to be zero. That is we take 

oj r (r, t = 0) = G(r, t = 0)=0. (24) 

Under these conditions we find 
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and B (r) is the baryon density, calculated from (||) with the initial pion field (22). 

Using the initial values we integrate the equations of motion to determine the fields at 
later times. As the fields propagate outwards, they diminish in size, making the non-linear and 
coupling terms in the Lagrangian less important. Thus we can define a radiation zone where the 
fields propagate as linear free massive fields. In the radiation zone, we can calculate the pion, 
f(k), rho, g(k), and omega, h(k), momentum distribution amplitudes from the expressions 
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respectively. Note that in extracting the fields we use that fact that they all satisfy p-wave 
equations in the asymptotic region. 

Let us examine the results of solving the dynamical equations. In Figure 1 we show the 
dimensionless pion field, F(r, t), as it propagates out to a time of 4 fm (c = 1). We see it begin 
as in (22) and diminish as 1/r. It does not disperse very much because the pion mass is small 
on the scale of energies of the problem. In Figure 2 we show a similar plot of the rho field. We 
plot the field, G(r,t)/gr reinstating the 1/gr factored out in (10). This field begins at t = as 
zero and builds from the non-linear couplings before beginning to diminish. The couplings and 
the large value of the p mass, make its behavior more complex than that of F. Figure 3 shows 
the radial part of the omega field, about which similar remarks apply. Figure 4 shows the time 
component of the omega field multiplied by r. This makes the plotted quantity dimensionless. 
We see that it is much smaller than F, and also that the factor of r makes it not decrease for 
large r, since the fields themselves all fall as 1/r for large r. 

A more intuitive picture of the time development of the classical fields comes from studying 



the evolution of the energy density. As we see from (21), that energy density has many parts. 
Also for our choice of initial condition all contributions to the energy density are spherically 
symmetric. Thus we do not plot the energy density as a function of r and t, but rather Airr 2 
times that density. The plotted quantity need then only be integrated with dr over all r to 
obtain the total energy. In Figure 5 we show the total radial energy density. We have checked 
that at each time slice the spatial integral of this density gives twice the nucleon mass. We see 
that the energy starts at t = rather compactly, propagates out and develops two "humps." 
The faster of these moves essentially along the light cone and is the energy in the prompt pions. 
The slower hump is the energy in the massive vector mesons and in the pions produced by the 
interaction of those vector mesons. Figure 6 shows the radial energy density in the pions only. 
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We see the prompt peak on the light cone but the slow peak is smaller than that in Figure 5 
since much of that slow peak is not in pions. Figures 7 and 8 show the radial energy densities 
for the rho and omega fields. These are both slower, because of the mass, and smaller (note 
the scale) than the corresponding pion energy. Figures 9, 10, 11 and 12 show various terms in 
the interaction energy. These are small (again note the scales) and, more importantly, decay 
quickly with increasing time or r. It is this diminished interaction energy at large r that allows 
us to define and asymptotic region where the fields propagate as free massive fields. It is these 
free fields that we use to construct the coherent state. 



5 Coherent State Calculation 

We use the asymptotic field strengths extracted in (E7|), ( p^ ) and ( p9| ) to construct quantum 
coherent states from the classical tt, p and to fields in the radiation zone. To obtain information 
on the pion multiplicity from these states, we need to project them into states of good four- 
momentum and isospin. We do this by known methods 16 1, that we have elaborated 



previously ||. First we construct the field operators(F, G and H) which create n, p and lu at 
the space time position x and pointing in the isospin direction T. They are 

F(x, T) = / d 3 k f{k)at ■ T e~ ik - x , (30) 



G(x, T) = I d A k g(k)bl ■ T e~ lk - x , (31) 



H{x) = / d 3 kh(k)cle- ik - x . (32) 



and 



Note that k ■ x = kot — k ■ x, with ko = y k 2 + M 2 ( M is the corresponding meson mass). 

From these we can form the coherent state with tt, p and u with 4-momentum and isospin 
projection as 

d 4 x dT 1 dT 2 . 
(2vr) 4 

xe iK-x^ e F(x,T 1 )+G(x,T 2 )+H(x) _ F ( X)Tl ) _ G{x,T 2 ) - H(x) - 1]|0> (33) 
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and we have defined a set of coupled harmonic functions of the ir and p isospin directions 

y IIz (T 1 ,T 2 )= ]T (II z \I 1 M 1 ;I 2 M 2 )Y IlMl (T 1 )Y hM2 (T 2 ). (34) 
liMr,l 2 M 2 



The states defined in ( p3| ) are orthogonal 

(I,I Z ,K\I' \I' Z ,K') = 5\K - K')6 n ,6 IzP l(K) (35) 
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where the normalization factor is given by 

[ d*x dT l dT 2 dT' l dT> 2 
1{K) ~ J (2^^ ^ynATi,T 2 )y Ih ( Tl ,T 2 ) 
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p g (x) = J d 3 p\g(p)\ 2 e-^ (38) 

and 

p h {x) = J d 3 p\h(p)\ 2 e-^ x . (39) 

It is easy to see that the normalization integral is independent of I z . We use the expansion 
method developed in ||, |6| to calculate the normalization integral 
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and 
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Similarly for the oj we have 
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The probability of finding pions, N p rhos and AT W omegas is given by 
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The probability of finding exactly n pions, from direct pions and from p and to decays, can be 
obtain from the joint probability ( |53|) 



P n = P{N*,N p ,N u ) (54) 

n=N w +2N p +3N u 

In the same spirit, we can also obtain the probability of having certain of number of mesons of 
each type and of given charge in the coherent state. That probability is 



p(N w+ ,N 7r -,N 7 ^, N p+ , N„- , Np , N w ) 



where 
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We now turn to the results from this formalism. 



6 Results 

One of the advantages of our approach to annihilation is that it contains very few parameters 
and that all annihilation channels are determined at once by those few parameters. Before we 
can discuss results, we must fix those parameters. The pion mass we fix at its experimental 
value of 140 Mev. The vector mesons, p and to, are assigned the same mass of m = 770 MeV. 
We choose the parameters, f n and g, to satisfy the KSFR relation |l4j m = \[2j^g. We choose 
three values for f v : 1) 93 MeV (experimental); 2) 75 MeV Q; 3) 62 MeV @- These last two 
values are popular choices in the Skyrmion literature where they are determined by getting the 
nucleon mass correct. In particular the middle value, f w = 75 MeV gives the observed nucleon 
mass in our formalism with gauged p and u fields |H| . The corresponding coupling constant g 
is: 1) 5.85; 2) 7.26; 3) 8.78. The only remaining parameter is the size of the initial spherical 
pion configuration. We fit that to obtain the correct single pion momentum distribution, giving 
a = 1.1 fm. This is certainly a reasonable range. All the remaining results to be discussed in 
this section follow from those few, fixed parameters. 

We begin with the single particle inclusive pion momentum spectrum from proton antiproton 
annihilation at rest. This is shown in Figure 13 where the data is compared with our calculation. 
As discussed above, one parameter in the initial pion configuration has been adjusted to obtain 
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the fit seen in Figure 13. Nevertheless we should stress that the entire machinery of our non- 
linear classical dynamics and subsequent quantum coherent state formalism intervenes between 
that initial configuration and the asymptotic momentum distribution seen in Fig. 13. That 
adjusting only one parameter in the initial configuration is capable of giving such a good fit 
to the distribution validates both the form of the initial configuration and the subsequent 
dynamics. In Figures 14 and 15 we show the corresponding momentum distribution for the rho 
and the omega. 

Nucleon antinucleon annihilation can occur in two isospin channels (I = 0, 1) with I z = 
for pp and I z = 1 for np. The mean number of mesons of each charge type for the three 
isospin channels are shown in Tables 1, 2, and 3, the different tables giving the result for the 
different choices of f n . The charged pions numbers shown are for the direct pions, while the 
final pion average number includes both direct pions and pions from vector meson decay. We 
see that the average number of all pions and the width of that distribution is rather insensitive 
to f n , but the details of the pion distribution and particularly the number of vector mesons 
is strongly dependent on f n . In all cases the average pion number and the width agrees well 
with the observed values. The results in Table 3 use the same parameters as in our previous 
calculation that did not include the p meson Q. We see that including the p has slightly 
reduced the number of direct pions, presumably to satisfy energy- momentum conservation. In 
the remainder of this section we use only the middle value of /„■ (/^ = 75 MeV) since it is 
the one fit to the nucleon mass in our classical field theory [pJ] ]. Thus apart from the range 
parameter, all our parameters are fixed by observed masses, pion, rho-omega and nucleon. 

In Figures 16 and 17 we show the pion number distribution for each of the isospin channels. 
These are both for primary pions and for those from vector meson decay. In each case we find 
the observed mean number of 5 with variance of 1. Also shown in those figures is a Gaussian 
distribution with the same average and variance. The Gaussian is nearly indistinguishable from 
the our calculation. Yet a Gaussian distribution is normally taken as evidence for a statistical 
or fireball process, while our process involves a rapid coherent pion wave. This should suggest 
a cautionary note to attempts to draw lessons about annihilation dynamics from the Gaussian 
pion number distribution. 

In Table 4 we show the branching ratios for the many resolved channels in proton antiproton 
annihilation at rest. Note that some of the channels involving more than two ir° are not resolved. 
Note also that the experimental signature of channels with vector mesons is uncertain and there 
are corresponding differences in the reported ratios from different experimental groups. The 
branches involving pions only in the second half of the table are for all pions, primary and 
secondary. We can easily convert our branching ratios for vector meson channels into pions 
since the vector mesons each have a principal pion decay mode; p ± — ► tt^ + 7r°, p° —* ir + + tt~ 
and u — > 7r + + 7r~ + 7T°. In Table 4 we show our calculation of the branching ratio for each of 
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the isospin channels separately, and then show the combined result based on an equal mixture 
of I = and 1 = 1. Our results would hardly change if we took the mix of 63% 1 = and 37% 
1 = 1 suggested by @. 

Also shown in the last row of Table 4 is the percentage of secondary pions. That is the 
fraction of pions coming from resonance decays. These are nearly all from p and u decay. The 
experimental fraction is 33%, and we find 30%, in quite reasonable agreement. The dynamical 
generation of vector mesons by our classical QCD is crucial to this result. 

The major features of the branching ratio data are correctly reproduced by our calculation 
with essentially no free parameters. The small channels come out small and the large large. 
For many of the large channels the agreement is excellent. Our biggest relative failing is for 
the 7r + 7r~ channel and for the 7r + 7r~7r channel. These are both two body channels since the 
7T + 7r~7r channel is dominated by Trp. We would expect the quantum corrections to be largest 
for the channels with the fewest quanta, and that is what we are seeing. Nevertheless the 
overall agreement lends considerable strength to the underlying treatment of classical QCD to 
generate the fields and a quantum coherent state projected on four momentum and isospin to 
give the various channels. It is the single quantum coherent state that is critical to explaining 
all the annihilation channels in a unified context. It is also important that four momentum 
conservation carry the constraints of phase space. These constraints are clearly seen in the 
broad-brush division into large and small channel. 

The effect of the classical field dynamics is most clearly seen in the generation of the observed 
vector mesons, since we start with zero vector meson fields. The first part of Table 4 shows 
that we are getting the vector meson channels qualitatively correctly and this is further verified 
in the correct fraction of secondary pions. Another way to see that our dynamics matters and 
that our agreement is not all from phase space is to do a pions only calculation. We construct a 
coherent state of pions only with an inclusive single particle momentum distribution fit to the 
empirically observed one |flQT| . We impose isospin and four momentum conservation and once 
again calculate the branching ratios to the pion channels. The results of this calculation are 
shown in Table 5. Note that for pions only there is an odd even isospin effect. That is states 
of odd numbers of pions can only contribute to I = 1 and even numbers to I = ||]. We see 
in Table 5 that both the mean number of pions and the details of the branching ratios are in 
far less good agreement with experiment than our previous calculation reported in Table 4 and 
including the full dynamics and the vector mesons. We conclude therefore, that just fitting the 
single pion spectrum and using it to construct a coherent state with phase space constraints 
imposed is not good enough. The underlying QCD dynamics also matters. 
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7 Summary and Conclusion 



We have shown that beginning with a classical field theory based in large Nc QCD we can give a 
good, nearly parameter free account of the major features of nucleon antinucleon annihilation at 
rest. Our premise is that annihilation occurs as a rapid process with a sudden burst of radiation 
and that that radiation can be described classically. We begin with the non-linear sigma model, 
the starting point of any large Nc QCD, and use gauge invariance to add classical rho and omega 
fields to the classical pion field. As these fields propagate out from the annihilation region, 
they decouple and can be well approximated by free fields. We use these asymptotic fields to 
construct a quantum coherent state and project that state onto states of good four momentum 
and isospin. We then calculate annihilation observables from that coherent state. We find good 
agreement with particle spectra and branching ratios. That agreement arises from our unified 
treatment of all the annihilation channels in terms of a single coherent state, which also has 
four momentum and isospin conservation imposed. It is interesting to observe that the number 
spectrum from the coherent states strongly resembles the experimental Gaussian distribution, 
while the conventional statistical approach (2(J (under the constraint of conservation of charge, 
isospin, etc.) failed to produce the Gaussian distribution. 

Our agreement with data is better than one might expect, given our simplified starting 
assumptions, but it is certainly not perfect. In particular the need for quantum corrections in 
the annihilation branches involving few final quanta is clear. The entire question of quantum 
corrections is an interesting and difficult one, but one that should now be seriously address 
in view of the demonstrated success of the first cut theory based on classical fields. We are 
presently studying quantum enhancements of the two pion channel to address the charge spin 



asymmetry seen in proton antiproton annihilation at higher energies [21, p2fl. Other leading 



quantum corrections are also under investigation. We are also investigating annihilation in 
flight. This would build in the two center nature of the process and allow a description of Bose 
Einstein correlations Further afield one might imagine a description involving classical 

SU (3) f fields to describe the small annihilation branches into K-mesons. 

The success in annihilation of the classical QCD starting point followed by quantum coherent 
states to describe the observed particles should give hope for similar methods in other problems 
where there are many pions but where the energies are low and nonperturbative QCD dominates 
the physics. Such approaches have already been developed for the disoriented chiral condensate 
H and have been suggested for heavy ion reactions All these situations raise the same 
questions. First is classical QCD appropriate for describing the dynamics and second can one 
estimate or calculate the leading quantum corrections? The difficulty of alternate approaches, 
the modest success reported here for annihilation, and the promised volume of high energy 
heavy ion data should make the investigation of these questions topical and fruitful. 
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Channel 


n,r+ n^o 


n p + n p - n p o 


n w 


n g 


I = I z = 


1.50 1.50 1.50 


0.09 0.09 0.09 


0.13 


5.44 0.81 


I=1I Z = 


0.99 0.99 1.95 


0.16 0.16 0.20 


0.13 


5.32 0.73 


7=14=1 


1.85 1.09 0.99 


0.29 0.06 0.16 


0.13 


5.32 0.73 



Table 1: Numbers of direct tt + , tt and 7r , p + , p and p ; a; for = 93 MeV. n is the total 
number of pions (direct and decay products) and a is the standard deviation calculated from 

Pn (0). 



Channel 


n n + n w - n^o 


n p+ n p- n p° 


n u 


n a 


1=04=0 


1.26 1.26 1.26 


0.06 0.06 0.06 


0.41 


5.33 0.86 


7=14=0 


0.80 0.80 2.09 


0.10 0.10 0.14 


0.33 


5.39 0.74 


1=14=1 


1.87 1.03 0.80 


0.20 0.04 0.10 


0.33 


5.39 0.74 



Table 2: Numbers of direct tt + , tt and it ; p + , p and p ; u for = 75 MeV. 



Channel 


n„.+ n,,- n^o 


n p + n p ~ n p o 


n u 


n a 


7 = I z = 


0.93 0.93 0.93 


0.02 0.02 0.02 


0.75 


5.18 0.82 


7=14=0 


0.61 0.61 2.24 


0.05 0.05 0.07 


0.55 


5.42 0.83 


7=14=1 


1.88 0.96 0.61 


0.10 0.02 0.05 


0.55 


5.42 0.83 



Table 3: Numbers of direct tt + , tt and tt ; p + , p and p , w for /V = 62 MeV. 
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Channel 


1 = 


Theory 
1 = 1 Combined 


Experiment 
CERN BNL 


TT°p° 


0.013 





0.006 


3.9 ±0.4 




0.013 





0.006 


1.9 ±0.3 







3.0 


1.5 


1.5 ±0.3 


5.8 ±0.3 


p°p° 


0.146 





0.073 


0.12 ±0.12 


0.4 ±0.3 


7T°U, 


0.0 


1.2 


0.6 


2.3 ±0.2 


0.7 ±0.2 


UJLO 


3.12 





1.56 


1.4 ±0.6 [ 


24 


1 


p°TT + TT-Tr° 


2.0 





1.0 


7.3 ± 1.7 


p^1T T TT + TT^ 


7.9 


10.3 


9.1 


6.4 ± 1.8 


LJTr + lT~ 


1.7 





0.8 


6.6 ±0.35 


3.8 ±0.4 


w2vr + 27r- 


4.3 





2.1 


1.3 ±0.3 


TT + TT~ 


0.02 


0.0 


0.01 


0.37 ±0.3 


0.32 ±0.04 




0.04 


0.6 


0.32 


6.9 ±0.35 


7.3 ±0.9 


27r + 2vr~ 


9.1 


3.0 


6.1 


6.9 ±0.6 


5.8 ±0.3 


27T + 27T-7r 


26.8 


19.8 


23.3 


19.6 ±0.7 


18.7 ±0.9 


37r + 37r _ 


13.8 


3.56 


8.7 


2.1 ±0.2 


1.9 ±0.2 


37T+37T-7r 


4.38 


0.61 


2.5 


1.9 ±0.2 


1.6 ±0.2 


n7T°, n > 1 


7.7 


15.7 


11.7 


4.1 ±0.4 


3.3 ±0.2 


7r + 7r _ n7r°, n > 1 


25.1 


39.8 


32.5 


35.8 ±0. 


8 


34.5 ± 1.2 


2vr + 27r-n7r°, n > 1 


12.8 


17.4 


15.2 


20.8 ±0.7 


21.3 ±1.1 


37r + 37r~n7T , n > 1 


0.03 


0.014 


0.022 


0.3 ±0.1 


0.3 ±0.1 


% of secondary its 


29.2 


31.3 


30.3 


33 



Table 4: Branching ratios, in percent, for proton antiproton annihilation at rest. Our calcula- 
tions are compared with experiments from [12]. We show each total isospin channel calculated 
separately. The "combined" column corresponds to equal mixture of I = and 1 = 1. In the 
last row we list the percentage of pions from the decay of rho and omega mesons. 
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Channel 


Pion-only Calculation 
1 = 1 = 1 Combined 


Experiment 
CERN BNL 




0.0004 0.0 0.0002 


0.37 ±0.3 0.32 ±0.04 


7r + 7r~7r° 


0.05 0.03 


6.9 ±0.35 7.3 ±0.9 


27r + 2vr- 


1.93 0.97 


6.9 ±0.6 5.8 ±0.3 


27r+27r-vr 


7.0 3.5 


19.6 ±0.7 18.7 ±0.9 


37T + 37T~ 


37.2 18.6 


2.1 ±0.2 1.9 ±0.2 


3vr + 3vr-7r 


10.5 5.2 


1.9 ±0.2 1.6 ±0.2 


n7r°, n > 1 


14 36 25 


4.1 ±0.4 3.3 ±0.2 


7T + 7r _ n7r°, n > 1 


16.8 30 23 


35.8 ±0.8 34.5 ±1.2 


27T + 27r~n7r°, n > 1 


20.9 16 18.5 


20.8 ±0.7 21.3 ±1.1 


37r + 37r _ n7r°, n > 1 


3 1.2 2.1 


0.3 ±0.1 0.3 ±0.1 


n, o" 


6.22, 0.83 6.4, 0.95 


5, 1 



Table 5: Branching ratios, in percent, for proton antiproton annihilation at rest, calculated from 
pions only, with experimental pion momentum spectrum. Data are compared with experiments 
from ||l2|| . We show each total isospin channel calculated separately. The "combined" column 
corresponds to equal mixture of I = and 1 = 1. In the last row we list the average and 
standard variation in the number of pions . 
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Figure 1: Pion field configuration F as a function of r and t. 

Figure 2: p = G/gr (in units of vector meson mass m) field as a function of r and t. 

Figure 3: uj r (in units of vector meson mass m) field as a function of r and t. 

Figure 4: u = tujq field as a function of r and t. 

Figure 5: Total energy density multiplied by Airr 2 (in units of MeV/fm) as a function of r and 
t. 

Figure 6: Pion energy density multiplied by 47rr 2 (in units of MeV/fm) as a function of r and 
t. 

Figure 7: Energy density associated with p multiplied by Airr 2 (in units of MeV/fm) as a 
function of r and t. 

Figure 8: Energy density associated with uj t multiplied by 47rr 2 (in units of MeV/fm) as a 
function of r and t. 

Figure 9: Energy density associated with p self-interaction multiplied by 47rr 2 (in units of 
MeV/fm) as a function of r and t. 

Figure 10: Energy density associated with 7rp interaction multiplied by 47rr 2 (in units of 
MeV/fm) as a function of r and t. 

Figure 11: Energy density associated with ituj interaction multiplied by 4-7rr 2 (in units of 
MeV/fm) as a function of r and t. 

Figure 12: Energy density associated with irpuj interaction multiplied by 47rr 2 (in units of 
MeV/fm) as a function of r and t. 
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Figure 13: Pion momentum distribution. The horizontal axis is the pion momentum. The 
vertical axis is the pion momentum distribution function, A-Kk 2 \f(k)\ 2 , see (|27|). The area under 
the curve gives the total direct pion number N^. The data is from fll2|. 



Figure 14: Rho momentum distribution. The horizontal axis is the rho momentum. The vertical 
axis is the rho momentum distribution function, A-Kk 2 \g(k)\ 2 , see (|28|). The area under the curve 
gives the total omega number N p . 



Figure 15: Omega momentum distribution. The horizontal axis is the omega momentum. The 
vertical axis is the rho momentum distribution function, 4irk 2 \h(k)\ 2 , see (|29|). The area under 
the curve gives the total omega number N w . 



Figure 16: Pion number distribution in / = channel. Solid squares are given by the distribu- 
tion P n (|53"|). Circles are the gaussian distribution with the same mean and variance. 



Figure 17: Pion number distribution in / = 1 channel. Solid squares are given by the distribu- 
tion P n (|53[). Circles are the gaussian distribution with the same mean and variance. 
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